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ABSTRACT 

Context. Recently measured straylight point spread functions (PSFs) in Hinode/SOT make granulation contrast in observed data and 
synthetic magnetohydrodynamic (MHD) data consistent. Data from earthbound telescopes also need accurate correction for straylight 
and fixed optical aberrations. 

Aims. We aim to develop a method for measuring straylight in the post-focus imaging optics of the Swedish 1-m Solar Telescope 
(SST). 

Methods. We removed any influence from atmospheric turbulence and scattering by using an artificial target. We measured integrated 
straylight from three diff'erent sources in the same data: ghost images caused by reflections in the near-detector optics, PSFs corre- 
sponding to wavefront aberrations in the optics by using phase diversity, and extended scattering PSF wings of unknown origin by 
fitting to a number of diff'erent kernels. We performed the analysis separately in the red beam and the blue beam. 
Results. Wavefront aberrations, which possibly originate in the bimorph mirror of the adaptive optics, are responsible for a 
wavelength-dependent straylight of 20-30% of the intensity in the form of PSFs with 90% of the energy contained within a ra- 
dius of 0.'6. There are ghost images that contribute at the most a few percent of straylight. The fraction of other sources of scattered 
light from the post-focus instrumentation of the SST is only ~10~^ of the recorded intensity. This contribution has wide wings with a 
FWHM -16'' in the blue and -34'' in the red. 

Conclusions. The present method seems to work well for separately estimating wavefront aberrations and the scattering kernel shape 
and fraction. Ghost images can be expected to remain at the same level for solar observations. The high-order wavefront aberrations 
possibly caused by the AO bimorph mirror dominate the measured straylight but are likely to change when imaging the Sun. We can 
therefore make no firm statements about the origin of straylight in SST data, but strongly suspect wavefront aberrations to be the 
dominant source. 

Key words. Instrumentation: miscellaneous - Methods: observational - Techniques: image processing - Techniques: photometric - 
Telescopes 



1. Introduction 

For many years, there has been a discrepancy between the con- 
trast in observed solar images and the corresponding images pro- 
duced by magnetohydrodynamic (MHD) simulations. It was not 
clear whether compensation for straylight in the observations 
or missing physics in the simulations were at fault. This has 
now been resolved by measurements of the point spread func- 
tion (PSF) in the Solar Optical Tele scope (SOT) on the Hinode 
spacecraft ( Wedemeyer-Bohm'2008 ) and compa rison with arti - 
ficial data (Wedemeyer-Bohm & Rouppe va n der Voort| |2009l). 
|S charmer et al.| ( |201Q[ see their introduction for a full account of 
the problem) recently initiated a project to measure the straylight 
sources of the Swedish 1 -meter Solar Telescope (SST; Scharmer 
et al.|2003a| to also correct SST images for the missing contrast. 



Scharmer et al.| (|20T0^ found that a significant part of the 



contrast reduction can be explained by high-order modes in the 
pupil wavefr ont phase that are not corrected by the Adaptive 
Optics (AO; Scharmer et al.|2003bj ) or image restorations with 
multi-frame blind deconvolution (MFBD; |Lofdahl||2QQ2] ) tech- 
niques because of the finite number of modes used in those tech- 
niques. The authors also implemented a method for correcting 
the contrast by means of the known statistics of atmospheric 
turbulence and simultaneous measurements of Fried' s parame- 



ter ro from a wide-field wavefront sensor. However, the contrast 
correction they found is not sufl&cient to reach the contrasts ex- 
pected from MHD simulations and Hinode observations. 

In this paper we continue the search for the missing con- 
trast by examining the post-focus optics of the SST. By inserting 
targets in the Schupmann focus, which is located close to the 
exit of the vacuum system, we removed all eflfects upstream of 
that point (atmosphere and telescope) and considered only error 
sources on the optical table. We estimated the wavefront PSF by 
using phase diversity (PD; Gonsalve s & Chidlaw|1979 Paxman 
fet al. 1992 ; Lofdahl & Scharmer 1994| ) and fitted the remaining 
straylight to a scattering PSF. In particular, we were interested 
in the integrated fraction of scattered light in the image data and 
the width of the scattering kernel. 



2. Imaging model 

We model the image formation as 



g = f^s^^SK + d + n 



(1) 



where g is a flat-fielded and dark-corrected data frame, / is the 
object, ^0 is the point- spread function (PSF) of wavefront er- 
rors 0, sk is a scattering PSF with extended wings, d represents 
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Table 1 : Scattering kernels 



Kernel name Definition 



FWHM 



Gauss 

Lorentz 

Moff^at 



KQ{r, cr) oc exp(-rV(2cr2)) 

/rL(r;r)oc(l+rVrV 
KM{r\a,p)oc{\+r^la^)-P 



= cr • 2 V21n2 



Wm = or • 2 V2i/^ - 1 



Voigt Ky{r, cr, 7) oc KG{r\ a) * Ki^{r\ y) Wy ^ 0.5346Wl + -^0.2166Wl + 

Notes. All kernels are functions of the radial coordinate r = (x^ + y^y^^. They are n ormalized in the Fourier domain by division with the value in 
the origin, is from|Mofl^(|1969j). The approximate expression for Wy is given by |01ivero & Longbothum| { [l977| ) and is claimed to be accurate 
to within a few %c. 
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Fig. 1: Straylight target drawing. Positions (x,y) in mm, diame- 
ters (0) in jim with tolerances of 1-2 jim. 



a residual dark level that was not subtracted properly in flat- 
fielding, and n is additive Gaussian white noise. These quanti- 
ties are all functions of the spatial coordinates (x,y), suppressed 
when possible for compact notation. The symbol ^ denotes con- 
volution. 

The scattering PSF, sk, is modeled as 



Sk = cS -\- (I - c) ' K 



(2) 



where S is the Dirac delta function, ^ is a convolution kernel 
equal to one of Km, Ki^, Kq, or Ky as given in Table [T] Because 
it is part of the scattering PSF, we will refer to ^ as a scattering 
kernel. Both the S function and the kernels are normalized in the 
Fourier domain by division with the value in the origin. Owing 
to this normalization, (1 - c) expresses the integrated fraction of 
this type of straylight. 



3. Data 

The data used in this experiment were collected with the SST 
( [Scharmer e"raL]|2003a| ) on 29 May 2010 between 17:10 and 
18:40 UT. 

The primary optical system of the SST is a singlet lens with a 
focal length of 20.3 m at 460 nm and a 98-cm aperture. A mirror 
at the primary focus reflects the light to a Schupmann corrector, 
which forms an achromatic focus next to the primary focus. In 
this focal plane we placed an artificial target to isolate the opti- 
cal aberrations and scattering downstream from this point. The 
target, manufactured by Molenaar Optics, is a 25 jim thick metal 
foil with six holes of diflTerent sizes, see Fig. [T] The manufac- 
turing tolerances are very small, but close inspection of the im- 
ages revealed small irregularities at the edges of the larger holes, 
probably caused by dust particles. Using a motor stage, we were 
also able insert a pinhole array used mainly for alignment pur- 
poses. 

The setup following the Schupmann focus is illustrated in 
Fig.|2] The beam expands (via the tip-tilt mirror) to a pupil plane 
at the location of the bimorph deformable mirror (DM). Here, the 



DM 




Fig. 2: Setup schematics. Light from the telescope enters through 
the target from upper left. TM = tip-tilt mirror; DM = de- 
formable mirror; RL = reimaging lens; DC = dichrocic beam- 
splitter; CT = correlation tracker; WFS = wavefront sensor; LCs 
= Liquid Crystals; FPI = Fabry-Perot interferometer. Only the 
cameras that produced images for this investigation are labeled 
(WB, PD). The figure is not to scale and the following inten- 
tional errors are introduced to save space: the angle to the beam 
reflected oflT the DC is drawn too wide and the transmitted and 
reflected beams from the first red beamsplitter are switched. 



telescope pupil is re-imaged by a field lens located just in front 
of the Schupmann focus. The pupil diameter is 34 mm at this 
location. There is a 35 mm diameter pupil stop at the DM. This 
defines the pupil for pinhole images. The reimaging lens then 
makes a F/46 beam parallel to the optical table. The light is then 
split by the 500 nm dichroic beamsplitter into a blue beam and 
a red beam. Both beams have several cameras behind diflTerent 
filters. Table [2] summarizes the cameras and setup parameters 
used for this experiment. 

In the blue beam one camera was nominally a "focus" cam- 
era and the other a defocused "diversity" camera of a PD pair. 
These two cameras and their beamsplitter were mounted on a 
common holder and could be moved together along the optical 
axis. This made it possible to generate additional focus diver- 
sities without changing the relative diversity between the wide- 
band (WB) and PD cameras. The holder and its cameras are cov- 
ered by a box that blocks light from directions other than the 
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Table 2: Setup in the two beams. 



Item 



Blue beam 



Red beam 



Wavelength (nm) 

No. of cameras 

No. of cameras used 

Cameras 

FOV (pixels) 

Image scale (arcsec/pix) 

Exposure time (ms) 



395.4 

4 

2 

MegaPlus II es4020 

2048x2048 

0.034 

10 



630.2 

3 

1 

SarnofF CAMIMIOO 

1024x1024 

0.059 

17 



beam. The relative rotation of the field in the WB and PD cam- 
eras was measured by comparing the grid patterns of the pinhole 
array images. The rotational misalignment is smaller than 0?1. 

In the red beam we used a single camera that was also cov- 
ered by a box to block external straylight. Here, we could only 
add diversities by moving the camera. 

The AO was running in closed loop on the central pinhole 
(D) of the straylight target while the artificial target data were 
collected, so aberrations are assumed to be stable during the 
data collection!^ The telescope pointing was moving near disk 
center to average out the granulation structure as well as pos- 
sible. Many exposures were collected at every position with all 
cameras, 3500 exposures in the blue and 10000 in the red, for 
the target data as well as for dark frames and flat fields. The 
high number of frames was necessary because we needed a high 
signal-to-noise ratio (SNR) to measure the very faint straylight 
far away from the pinholes. The target data were corrected for 
bias and gain variations. The corrected image gi is calculated as 



g. = ih,-d)l{f-d\ 



(3) 



where hi is the observed data frame, d is the average dark frame, 
and / is the average flat field image. To minimize the influence 
from solar features in the averaged flat-field images, the individ- 
ual frames were collected while the telescope was circled over a 
quiet area near disk center and the DM produced random wave- 
fronts. 

The images collected in diff'erent focus positions were not 
recorded simultaneously and we took steps to ensure that tiny 
image movements did not violate the assumption that we made 
in PD processing, i.e. that there is a common object in all im- 
ages after summation. In a first pass, the data frames gi were 
co-added in their original form. In a second pass, each exposure 
was aligned to the summed WB image at nominal focus (WBO 
image) from the first pass (separately in blue and red) with sub- 
pixel precision through centroiding (center of mass) on one of 
the holes, and a new sum was formed. By using the same image 
as an alignment target, both the individual images in each sum 
and also the sums at diff'erent focus positions were well aligned. 

In the remainder of this paper we will refer to these summed 
images simply as the images, corresponding to g in Eq. ([T]). The 
WBO images are shown in Fig. [3] 



^ Because the target was repeatedly removed from the beam (when 
collecting flat fields) and inserted (when collecting target data), the mir- 
ror was alternately heated by the sunlight and allowed to cool again. 
Because these thermal variations cause changes in the mirror's shape 
and stress, so the level to which the assumption is valid is difficult to 
predict. 




Fig. 3: Straylight target images, WBO in both beams, displayed 
with the same log scale. Left: blue beam; right: red beam. 
Compare Fig.[T] 

4. Straylight measurements 

4.1. Ghost images 

Figure [3] shows a variety of weak duplicates of the pinhole im- 
ages, commonly referred to as ghost images. The most signif- 
icant ghost images are shifted by only a few arcseconds and 
appear to be at diff'erent focus positions. These contributions 
must come from reflections in the beam splitters and other optics 
within a few cm from the detector focus, and should therefore be 
present during ordinary observations of the Sun as well. There 
are many weaker contributions that appear to be mirror images 
in the blue and repetitions in the difl'erent taps of the red CCD, 
probably originating in the camera electronics]^ 

The ghost images are not included in the image formation 
model of Eq. ([T]), but we can measure them here and then mask 
them in later processing. The strongest more or less in-focus 
ghosts are ^ 1% in both the blue and red (measured on hole 
A). 

4.2. Wavefront aberrations 

We estimated the wavefront error PSF, s^p, through PD with the 
MOMFBD C++ code of |van Noort et al.| ( |2QQ 5). We processed a 
subfield centered on the 20 jim straylight hole F. In the blue we 
used 256x256 pixels and in the red 128x128 pixels, the diff'erent 
sizes are necessary because of the diff'erence in image scale. The 
top row of Fig. |4] shows the central part of this subfield. We can 
see directly that the chosen focus positions were not ideal. PDO 
is very similar to WB- and PD+ to WBO in the blue, so instead 
of six diff'erent diversity channels, we have in reality only four. 
The focus is somewhere between WBO and WB+ in the red; it 
is probably preferabl^to have one channel close to the focal 
plane, to constrain the object estimation. 

Careful pre-processing is easier with point-like targets than 
with extended scenes, which facilitates inversions. In particular, 
it is possible to fine-tune the intensity bias (dark level) correc- 
tion, and then make sure the total energy is the same in all focus 
positions. Such refined dark-level correction was performed for 
the purpose of PD processing. First, all images in a PD set were 
normalized with respect to the average intensity in the inner part 
of the 1 mm target hole (A). We assumed this to be the best avail- 
able measure of the intensity level, insensitive to the blurring of 



2 The Sarnoff CAMIMIOO CCD is organized in 2 by 8 taps (sub- 
fields), that are read out in parallel. 

^ We are not aware of any studies about optimal distribution of sev- 
eral focus diversities. 
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Fig. 4: PD results, the central parts of target hole F images in log scale. Top row: observed image; center row: recreated images 
M = 231; bottom row: recreated images M = 36. The FOV shown is 2!2x2!2 in the blue and V!9xV!9 in the red. Beam, camera 
and diversity as noted in subcaptions a)-i). 



the edge of the hole. The WBO images, which were collected 
in nominal focus, we then subtracted a dark level measured as 
the peak (as defined by a Gaussian fit) of the histogram of im- 
age intensity within the PD processing subfield. This corrects for 
some of the straylight from hole A. For the data from other focus 
positions, we subtracted a dark level corresponding to the diff'er- 
ence in median intensity over the entire frame (defining the dark 
level). Finally, the images were normalized to the same mean 
value, making the total energy the same in all focus positions. 
The dark corrections, d, in this refinement step were on the or- 
der ~10~^ of the intensity in hole A in the blue and almost ~10"^ 
in the red. 

In the PD processing the wavefronts were expanded in at- 
mospheric Karhunen-Loeve (KL) fun ctions, express ed as linear 
combinations of Zernike polynomials fRodd ier|1990| ). These KL 
functions are ordered as the dominating Zernike polynomial in 
the notation of Noll (1976 ), and not by monotonically decreas- 
ing atmospheric variance. Although we were not measuring at- 
mospheric wavefronts, we chose to use KL functions because 
bimorph mirrors naturally produce KL modes. We assumed that 
the alignment of the PD channels described in Sect. |3] is suf- 
ficient for PD processing, and accordingly did not include tilt 
coefficients in the fit. We defined M as the index of the highest- 
order mode used, i.e., we used KL modes 4-M. 

In addition to the KL parametrization of the unknown wave- 
front, we also estimated the (Zernike) focus diversities with re- 
spect to the WBO images, starting from the nominal values. 
Because the WB and PD cameras sit on a common mount and 
are therefore moved together, a single additional focus shift had 
to be determined when the three images from the PD camera 
were added. The diversities estimated with diff'erent M vary by 
no more than a few tenths of a mm and we used the values es- 
timated with M = 210 for both M used here. The nominal and 
estimated diversities are shown in Table |3l 

The magnitudes of the estimated diversities are all lower than 
their nominal values. This could be caused by a mismatch in 
the image formation model. We used the diameter of the ordi- 
nary telescope pupil, re-imaged on the bimorph mirror. In real- 
ity, there is diff'raction in the pinhole, making the re-imaged pupil 
slightly fuzzy. The 7 mm defocus of the PD camera is particu- 
larly underestimated. This may be because of a mistake in the 



Table 3: Focus diversities. 



Diversity 




Blue 




Red 


WB- 


WB+ 


PDO 


WB- WB+ 


Nominal (mm) 
Estimated (mm) 


-4.0 
-3.8 


+4.0 
+3.7 


+7.0 
+4.8 


-7.0 +7.0 
-6.2 +5.4 



Notes. Diversities are expressed in mm shift along the optical axis, rel- 
ative to WBO in each beam. The PD+ diversities in the blue are simply 
the sums of the WB+ and PDO diversities. 

Table 4: RMS and Strehl ratio of the wavefront phase. 



Beam 


M 


Measured 




In focus 




(waves) 


(waves) 


^^obs 


Cf 500 

(waves) 


^500 


Blue 


36 


0.059 


0.059 


0.87 


0.047 


0.92 


Blue 


231 


0.105 


0.103 


0.66 


0.082 


0.77 


Red 


36 


0.100 


0.052 


0.90 


0.065 


0.85 


Red 


231 


0.121 


0.070 


0.82 


0.088 


0.74 



Notes, cr^^^^ is the RMS of the PD-estimated wavefront in waves at 
the observations wavelengths, given for the entire estimated wavefront 
phase, as well as after subtracting the best-fit Zernike focus component. 
Strehl ratio calculated as = exp{-(2;rcr^)^). Subscript 500 denotes 
A = 500 nm. 



setup, so this camera was in fact defocused by less than 7 mm. 
These estimates confirm that the diversity setting of the PD cam- 
era was not very diff'erent from the step size of the (-, 0, +) po- 
sitions. 

Along with the observed images, g, in Fig. |4]we also show 
the recreated images, g = f^S(p, based on the PD estimates / and 
^, for diff'erent M. An inspection of these images clearly shows 
that the PSFs must be satisfactorily estimated, particularly in the 
blue M = 231 case. 

The estimated wavefront phases without focus are shown in 
Fig. [5] The M = 231 estimates in the blue and red show simi- 
larities, although the red wavefronts are less well resolved than 
the blue wavefronts because of the smaller subfields used and 
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(a) Blue 36 



(b) Blue 231 



(c) Red 36 



(d) Red 231 



Fig. 5: Estimated wavefronts, 0, in focus. Beam and M as noted 
in subcaptions a)-d). The scaling is between the same min and 
max for the two wavefronts from the same beam. 



possibly because of the smaller number of images with different 
diversities. They both have two bright rings and a dark gradient 
at perimeter. Both rings have brighter bumps at approximately 
the same positions. Particularly in the blue wavefront, the pat- 
tern resembles that of the electrodes on the bimorph mirror. The 
M = 36 wavefront estimates are much less resolved. Some sim- 
ilarities to the 231 -estimates can be seen along the outer bright 
ring. However, with the smaller M, the KL functions apparently 
cannot represent the inner ring pattern visible in the M = 231 
wavefront estimates. 

The RMS values of the estimated wavefronts are listed in 
Table [4] The WB camera in the "0" position was out of focus by 
approximately 0. 1 rad in the blue and 0.5 rad in the red. To make 
the RMS comparable between beams, we also show the RMS 
after subtracting the best-fit focus contribution as well as cal- 
culated for a common wavelength, 500 nm. This makes the 231- 
estimates in red and blue agree, (T500 = 0.085+0.003 waves. The 
36 and 231 -estimates do not agree in either beam. We are more 
confident in the blue 231 -results because 1) the results agree be- 
tween blue and red, 2) the details in the wavefront are the so- 
lution converged to when M is increased in increments from 36 
to 231, and 3) the solution give estimated PSFs that recreate the 
observed data well, see particularly the large-diversity images in 
Fig.|4] The Strehl ratio for this solution (see Table]?]) is 0.66 in 
the blue. Taking the two 231 solutions in red and blue together 
we obtain a Strehl rati(f]of 75 ± 2% at 500 nm after removing 
the focus error. 

Instrumental contributions to the wavefronts can in princi- 
ple be estimated and corrected in solar data with the normal 
MOMFBD image restoration, except that we usually use 36 
modes, far less than 231. The instrumental high-order modes are 
not included in the ro measurements of the wide-field wavefront 
sensor s o they are also no t dealt with in the high-order compen- 
sation of [Scharmer et al. ([2010). They are therefore an indepen- 
dent contribution to the lowered contrast in our solar images. 

In Fig.|6]we show the enclosed PSF energy as a function of 
the radial coordinate. The enclosed energies for ^0 and reach 
90% at r = a.' 17 and a.' 62, respectively, in the blue and at a.' 27 
and 0^.^55 in the red. 



4.3. Extended wings 

We attempted to fit a scattering PSF with extended wings to the 
WBO images shown in Fig. [3] We preferred to use the 1 mm 
hole A, because it transmits more light than the other holes and 
therefore has better SNR in the far wings. To remove the pollu- 
tion from the other holes and from the ghost images, we defined 
binary masks to be used when fitting, see Fig. [7] 

^ The Strehl ratio is the observed peak of the PSF divided by the peak 
of the theoretical PSF for a perfect imaging system. 
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Fig. 6: Enclosed energy for the diff'raction-limited PSFs, sq, and 
for the estimated wavefront aberration PSFs, S(p. The dotted lines 
illustrate the radii for 90% enclosed energy. 




Fig. 7: Masked stray light target images, displayed with the same 
log scale as in Fig. [3] Left: blue beam; right: red beam. 



For the WBO images used for fitting the straylight model, the 
second-pass sums were corrected for a pattern of stripes running 
in the y direction. This correction was made by subtracting from 
all rows a smoothed version of the mean of rows 1-10. The av- 
erage of this correction was -10"^ in the blue and almost -10"^ 
in the red, representing a correction of the dark level similar to 
the levels corrected in the previous section. (Fig. [3] shows the 
images after this step.) For this step the images were also nor- 
malized with respect to the average intensity in the inner part of 
the 1 mm hole (A). 

A binary representation of hole A was generated by thresh- 
olding the WBO image at 50% of the maximum intensity and 
removing contributions from the other holes. Artificial images 
were then made by convolving this binary image with the appro- 
priate PSFs. 

Figure [8] shows an azimuthal average (taking the mask into 
account) of hole A, zoomed in on the hole perimeter. (The zoom 
makes a discrepancy in image scale of about 1% apparent.) We 
show the artificial hole convolved with three diff'erent PSFs and 
the observed hole. The artificial image made with PD estimated 
match the observed data very well in the blue. Comparing 
with the diff'raction limited image, it is obvious that we need to 
model the wavefront PSF to isolate other sources of straylight. 
The match is fairly good also in the red but the near wings are 
overestimated. The two PD estimated PSFs (M = 36 and M = 
231, resp.) give almost indistinguishable results for this purpose. 
The straylight we will try to model with the scattering PSF, sk, is 
the component that shows up at r ^ 5V8 in the blue and r ^ 6V2 
in the red. 

We now fitted Eqs. ([T]) and ^ to the data, using the M = 
231 estimates of as fixed contributions from the wavefronts. 
We used the Levenberg-Marquardt algorithm as implemented in 
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r/1" r/1" 

Fig. 8: Observed and artificial images, binary 1 mm hole convolved with and as measured with two diff'erent M. Left: blue 
beam; right: red beam. 



the MPFIT package for IDL ( |More|1978[ |Markwardt|2009l ). In 

practice we fitted with two parameters ci and C2 substituted for 
c and (1 - c) in Eq. ^ to allow for normalization diff'erences in 
real and artificial data. We then report c = ci/(ci + C2) as the 
result of the fit. 

The results of the fits are shown in Fig. [9] and in Table [5] 
With any of the kernels we get c > 0.99 in both beams, which 
appears to be a very robust result by the agreement between the 
diff'erent fits. In the blue we get 0.2-0.4% scattered light and in 
the red 0.1-0.4%. 

The d parameter is returned with low values, on the order 
of 10"^. This indicates that the subtraction of the average top 
rows worked well as a dark correction refinement. Note that in 
the red, outside the signal-dominated radii, there is just noise, 
seemingly around zero. In the blue there seems to be a remaining 
uncorrected dark level that is not modeled by the d fit parame- 
ter. This could be caused by a geometrical asymmetry, e.g., by 
insuflftcient masking of the other holes or a dark-level gradient. 
It does not seem to ruin the fits. 

We have less confidence in the details of the red fits than in 
the blue fits because the artificial image based on overesti- 
mates the core of the PSF, see the right panel of Fig. [8] However, 
it is apparent from Fig. [9] that there is less straylight in the red 
than in the blue. Comparison of the observed red data (black in 
the red plots) with the observed blue data (gray) reveals a factor 
of ~3 diff'erence. This would favor the lower estimate of 0.1% in 
the red and make the Lorentz kernel fit less believable. 

In the blue the fit with a Gaussian kernel seems to fail at large 
radii. It also has the largest FWHM, ~20''. The Lorentzian ker- 
nel follows the signal dominated curve better than the other ker- 
nels but the Voigt and Moff'at kernels give similar results. These 
all have FWHM In the red the Lorentzian kernel dif- 

fers from the other kernels by giving a larger FWHM. While the 
Moff'at kernel /S is estimated with a very high value, which is 
compensated for by an a that is also very large, and the Voigt 
kernel is almost degenerated to a Gaussian, the latter three ker- 
nels agree on a FWHM of 34'' . We conclude that the Voigt and 
Moff'at kernels, by virtue of their two fit parameters, are better 
able to represent the true shape of the extended wings. Using 
their results, we estimate the FWHM to W in the blue and 34'' 
in the red. 

There are noticeable residuals at 10''+ a few arcsec in the fits 
shown in Fig. |9] Particularly in the red, where this feature cor- 



responds to 36% of the total intensity! Because of the log scale, 
this residual looks much lower in the blue but the corresponding 
error in the blue is in fact 15%. This may sound alarming but it 
appears this is needed to cancel fitting errors within the hole and 
at the hole perimeter. There are several possible reasons for the 
bad fits at small radii. There are small imperfections in the shape 
of the hole that could cause diff'raction. The 50% threshold does 
not necessarily produce a binary representation of the hole with a 
correct size. The thickness of the metal foil may cause reflections 
in the interior walls of the hole. There are remaining variations 
in intensity of unknown origin within the hole. Paricularly in the 
red, the PD estimated s^p does not fit the core of the PSF. Owing 
to the truncation of the wavefront expansion, the estimated 
may also under-represent the near wings of the PSFs. Thermal 
relaxation may cause variations in the S(p errors around the hole 
perimeter during the data collection. Despite the limited accu- 
racy at small radii, the residuals for the extended wings at larger 
radii are low and we believe the kernel fits can be trusted because 
the 2D fitting is dominated by the larger areas where the radii are 
large. 

5. Discussion 

The wavefront aberrations are quite significant, corresponding to 
a Strehl ratio of -75%. The structure of the estimated wavefronts 
(M = 231) with their 6- and 12-fold symmetries suggests that the 
bimorph mirror of the SST AO is responsible. Wavefront aber- 
rations originating in the instrumentation will be at least partly 
corrected for by our standard MOMFBD image restoration. Part 
of the wavefront aberrations came from modes of higher order 
than we normally include in this processing and it is unclear how 
much of this would be corrected for. With the point-like object 
used here, the magnitude of the estimated wavefront strongly de- 
pends on the number of included modes, M. However, as shown 
by Scharmer et al. (2010), when the object is solar granulation, 
the MFBD/PD-type problem is less constrained and an estimated 
low-order wavefront tends to also represent the blurring caused 
by the higher-order modes that are not included. 

The aberrations give straylight that is mostly contained in 
the first diff'raction rings of 90% of the energy is within a 
radius of Off 6 (see Fig. |6]). The origin of this surprisingly high 
level of wavefront errors is unknown. One possible source is 
the mounting of the 1 mm thick deformable mirror, which is 
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Fig. 9: Fit results, angular averages of the masked images, centered on the 1 mm hole A. Beam and M as noted in the plot legends. 
Black: observed image; gray (in the red plots): blue observed image; brown: binary image convolved with estimated with M = 
231; Other colors: binary image convolved with as well as Sk fitted with kernels as listed. 

Table 5: Scattering kernel fit results 



Blue 



Red 







c 


Kernel parameters 


FWHM 




c 


Kernel parameters 




FWHM 


x' 


36 


Gauss 


99.8 


0-= 8.31 




19^.^6 


0.0143 


99.9 


cr- 13.74 






3274 


0.0595 


36 


Lorentz 


99.6 


7= 6.61 




1372 


0.0142 


99.7 


7= 17.78 






3576 


0.0595 


36 


Moffat 


99.7 


Qf= 9.03 


1.36 


1477 


0.0142 


99.9 


a = 416.03 




458.67 


32r4 


0.0595 


36 


Voigt 


99.7 


cr= 3.98 


7 = 4.23 


1477 


0.0141 


99.9 


cr= 13.64 


r = 


0.0008 


32ri 


0.0595 


231 


Gauss 


99.8 


cr= 8.77 




2a.'6 


0.0148 


99.9 


cr= 14.61 






34r4 


0.0593 


231 


Lorentz 


99.6 


7= 7.38 




14r8 


0.0147 


99.6 


7= 19.83 






39'.'7 


0.0593 


231 


Moffat 


99.7 


a = 11.04 


1.56 


i6r5 


0.0147 


99.9 


a = 524.28 




644.53 


34r4 


0.0593 


231 


Voigt 


99.7 


cr= 4.59 


7 = 4.26 


i6ri 


0.0147 


99.9 


cr= 14.39 


7 = 


0.0008 


3379 


0.0593 



Notes. M is the max KL index used in PD. The parameters a, 7, and cr are in units of arcsec, while /3 is dimensionless. The c parameter is the 
percentage of the intensity that is not contained within the far wings. 



clamped between two 0-rings with the tension adjusted manu- 
ally by 12 screws. Another possibility is high-order aberrations 
induced by a large focus error imposed on the mirror. 

pom see the supporting online material) 



Scharmer et al. 



compared umbra intensity and granulation contrast in 630 nm 
SST data to data from SOT/Hinode and inferred a 50% stray- 
light level with a small FWHM of less than T\ consistent with 
this stray light originating from aberrations. However, the aber- 
rations measured here can only account for about 1/3 of the 50% 



straylight and would have to be stronger by a factor 2 than those 
measured here to fully explain the contrast in the observed gran- 
ulation data. 

Our best estimate of the scattering PSF, sk, is that it gen- 
erates only 0.3% straylight in the blue (-397 nm) and 0.1% in 
the red (~630 nm). In our best fits, the FWHM of the scatter- 
ing kernel K is ~16'' in the blue and ~34'' in the red, using the 
Voigt and Moffat kernels. The Gauss and Lorentz kernel fits give 
similar results. 
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It is common to truncate the wings of the scattering kernels 
at some radius but our fits work well without truncation, possibly 
because we are including the dark level d in the fit. We cannot 
exclude that some pollution from the holes B-F influenced the 
sk fits. A better stray light target should have only the 1 mm hole 
and the smallest hole. 

After the ordinary dark-current calibrations, which consist of 
a subtraction of an average dark frame, the dark level was ma- 
nipulated in various ways. For PD processing, we estimated and 
subtracted ~10"^ of the intensity in hole A from the subfield con- 
taining hole F, measured basically as the modal value within the 
PD processing subfield. Before the straylight kernel fitting, we 
also subtracted on the order of -10""^ based on the top 10 rows of 
pixels. Because the kernel fitting involved a dark level, J, which 
was estimated to -10"^, we can be reasonably sure that the stan- 
dard dark-correction is correct to the 10"^ level. However, this 
is insignificant in comparison with the ghost images, the most 
significant of which could be measured to ~1% of the hole A 
intensity. Consequently, we have to allow for an additive com- 
ponent at the -10"^ level, which is not modeled by the scattering 
kernels. 

6. Conclusions 

We have proposed a procedure for measuring the amount of 
straylight in the SST post- focus instrumentation and applied it 
to one wavelength in the blue beam and one wavelength in the 
red beam. The strength of the method is that we simultaneously 
measured wavefront aberrations (with phase diversity methods) 
and "conventional" scattered light from the same data. Thus, we 
were able to separate these two sources of straylight. 

The dominant contributions to straylight are high-order aber- 
rations, causing a reduction of the Strehl ratio. The estimated 
Strehl ratio at 390, 500 and 630 nm is 66%, 75%, and 82%, 
resp., corresponding to integrated straylight within a radius of 
~0'.'6 on the order of 34%, 25%, and 18%, resp. The second- 
most important source of straylight is multiple weak, out-of- 
focus ghost images. The combined eff'ect of these is difficult to 
estimate but most likely contributes less than a few percent at all 
wavelengths. The smallest contribution measured is from "con- 
ventional" straylight in the form of a PSF with very wide wings 
(16'' in the blue, 34'' in the red), this is estimated to contribute 
0.3% in the blue beam and 0.1% in the red. 

The wavefront contribution could potentially be higher when 
observing the Sun because of the heat that aff'ects the mirror. 
Considering the small number of optical surfaces and the high 
optical quality of the SST telescope optics, other main contribu- 
tions to the lowered contrast most likely are of atmospheric ori- 
gin. This could come partly from uncorrected high-order wave- 
front aberrations, as discussed by Sch armer et al. (2010|, and 
partly from scattering by dust particles in the Earth's atmo- 
sphere. 
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of La Palma by the Institute for Solar Physics of the Royal Swedish Academy 
of Sciences in the Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 

Note added in proof In October 2011 we were upgrading the 
wavefront sensor of the SST AO system. We found that switch- 
ing off' the mirror voltages results in a change in the science 
focus by approximately 9 cm along the optical axis. Due to 
limited electrode resolution, deformable mirrors cannot produce 
even low-order modes perfectly. When such a large focus has to 
be compensated for, there are therefore unavoidable high-order 
wavefront errors. We have simulated our setup and calculated 
the wavefront phase corresponding to the diff'erence between a 
9 cm defocus and the approximate focus produced by the mirror 
in order to compensate for it. The resulting RMS wavefront is 
0.14 waves at 500 nm and the corresponding Strehl ratio is 0.48. 
This is even worse than the eff'ects reported in Table|4]by a factor 
1.6 in RMS wavefront. 

We note again that these measurements (and simulations) are 
not directly comparable to solar observations. Nevertheless, this 
eff'ect may well contribute by a signiffcant amount to the loss 
of contrast in SST data. The SST optical setup will be modiffed 
to take this into account when we install our new 85 -electrode 
deformable mirror during the summer or fall of 2012. Before 
that, we will also try to compensate for the eff'ect by introducing 
a corresponding focus change in the Schupmann system. 
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